PROCEEDINGS 

°F SCIENCE 



ANL-HEP-CP-13-4 



Endpoint Logarithms in e^e~ J / \i/ -\- r] ^ 



Geoffrey T. Bodwin* 

HEP Division, Argonne National Laboratory 
E-mail: |gtb@hep . anl . gov 



Hee SokChungt 

Department of Physics, Korea University 
E-mail: |neville@korea . ac . kr 



Jungil Lee 

Department of Physics, Korea University 



E-mail: jungil@korea . ac . kr 



We investigate the origin of the double logarithms of /ni^ that appear in the calculation of the 
cross section for e+e^ Tjc at next-to-leading order in the strong coupling a^. Here, 

is the square of the center-of-momentum energy, and is the charm-quark mass. We find that, 
diagram-by-diagram, the double logarithms are accounted for by Sudakov double logarithms and 
endpoint double logarithms. The Sudakov double logarithms cancel in the sum over all diagrams, 
but the endpoint double logarithms do not. We reinterpret the endpoint double logarithms in 
terms of a leading region of loop integration in which a spectator fermion line becomes soft and 
collinear This reinterpretation may simphfy the process of establishing an all-orders factorization 
theorem for this helicity-flip process, which, in turn, might allow one to resum logarithms of 
/m^. to all orders in a,. 



Xth Quark Confinement and the Hadron Spectrum 
8-12 October 2012 

TUM Campus Garching, Munich, Germany 



*Work supported by the U. S. Department of Energy, Division of High Energy Physics, under Contract No. DE- 
AC02-06CH 11357 
' Speaker. 



© Copyright owned by the authoi (s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Endpoint Logarithms in e^e J/\j/+ rjc 



Hee Sok Chung 



1. Introduction 



The exclusive double-chamionium production cross section (j{e^e^ — )• J/y+Vc), which has 
been measured by the Belle ^ and BABAR [^] collaborations, has provided serious challenges 
to the nonrelativistic QCD (NRQCD) factorization formalism ||, ^, ^, ||]. The theoretical value 
for the cross section at leading order (LO) in the strong coupling a, and the heavy-quark velocity 
V is almost an order of magnitude smaller than the measured rate. The discrepancy between 

theory and experiment seems to have been resolved by the inclusion of higher-order corrections, 
which include QCD corrections of next-to-leading order (NLO) in a, |11, 12] and corrections of 
higher order in v (relativistic corrections) [13, 14, 15 1. Here, v is the velocity of the heavy quark 
(Q) or heavy antiquark {Q) in the quarkonium rest frame. 

As the authors of Ref. JT^ ] have pointed out, the large NLO QCD coiTections involve double 
logarithms of /m^, where is the square of the e+e^-center-of-momentum energy and is the 
charm-quark mass. These logarithms are sufficiently large that it may be necessary to resum them 
to all orders in in order to obtain a reliable theoretical prediction. Typically, the resummation 
of large logarithms is carried out through a factorization of the contributions that arise from small 
momentum scales from the contributions that arise from large momentum scales. A first step in 
deriving such a factorization is to identify the regions of loop-momentum integrals that produce the 
large logarithms at the leading nonzero power of [p^]. 

In this paper, we identify the leading loop-momentum regions that give rise to the double log- 
arithms in the NLO QCD coiTections to the process e+e^ J/y + Tlc- According to the NRQCD 
factorization formalism [Q], the amplitude for this process at LO in v, can be written as a product 
of a short-distance coefficient with the NRQCD long-distance matrix elements (LDMEs) for the 
evolution of QQ pairs into the J/y and the rjc. The same short-distance coefficient appears in the 
amplitude for the production of two QQ pairs that have the same quantum numbers as the corre- 
sponding quarkonia. The short-distance coefficient can be obtained perturbatively by comparing 
the fuU-QCD amplitude i£/[e+e- QQi^Si) + QQi{^Sq)] with the NRQCD amplitude, which 
consists of the short-distance coefficient times the NRQCD LDMEs for the QQ states. (Here the 
subscripts 1 indicate that the QQ pairs are in color-singlet states.) The double logarithms of Q^/m^. 
arise solely from the full-QCD amplitude because the NRQCD LDMEs for the QQ states are in- 
sensitive to momentum scales of order m^. or larger. In Sec. ^, we carry out the calculation of the 
double logarithms of Q^ /ml at NLO in for the process e^e^ -^J/y+ by examining the NLO 
QCD coiTections to the full-QCD process. We find that the double logarithms arise from Sudakov 
or endpoint regions of loop momenta, and we identify the contribution that arises from each region 
for each NLO Feynman diagram. Our results for the double logarithms agree with those that were 
obtained in the complete NLO calculations |11, 12]. We find that the Sudakov double logarithms 
cancel in the sum over diagrams and that power-divergent contributions from the endpoint region 
vanish. In Sec. ^, we give a general analysis of the Sudakov double logarithms that elucidates the 
reason for their cancellation in the sum over diagrams. In Sec. ^ we give a general analysis of the 
endpoint region that establishes the absence of power-divergent contributions. We summarize our 
results in Sec. S. 
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2. Calculation of double logarithms 

In this section, we evaluate the double logarithms that appear in the NLO QCD corrections to 
the amplitude e^e^ J /y + Vc^ and we identify the momentum regions that are associated with 
the logarithms. The process e^e^ — )• QQi{^Si) + QQi CSq) is composed of e^e^ — 7*, followed 
by 7* — )• QQi {^Si ) + QQiCSo)- Because the process e^e^ — )• 7* does not receive QCD coixections 
in relative order a^a^, we need to consider only the process 7* — )• QQi{^Si) + QQi{^So) in the 
evaluation of the NLO QCD corrections. 

2.1 Kinematics, conventions, and nomenclature 

Now we describe the kinematics, conventions, and nomenclature that we use in calculating 
the double logarithms and throughout this paper. We work in the Feynman gauge. We use the 
light-cone momentum coordinates k = [k'^,k^,k_i] = [{k^ + k^)/\/2,{k^ — k^)/V2,kj_] and work 
in the e+e^-center-of-momentum frame. Because our calculation is at LO in v, we set the relative 
momentum of the Q and Q in each charmonium equal to zero. Then, the momenta of the Q and 
Qin the //i/A are both p= [{^/P^^Tn^ + P)/\^,{^/P^^Tn^ - P)/V2,0±], and the momenta of 
the Q and Q in the rjc are both p = [{y^W+n^ - P) / V2,{y^W+n^ + P) /V2,0_i], where P is 
the magnitude of the 3-momentum of any of the Qs or Qs. The momentum of the virtual photon is 
Q = 2{p + p), which implies that = 16(P^ + ml). If a momentum k has light-cone components 
whose orders of magnitude are PX[l, (t7+)^,T7+], then we say that k is soft if A ^ 1, and we say 
that k is collineai- to plus if T7 + <C 1 . If ^ has light-cone components whose orders of magnitude 
are PA[(t]^)^, 1, T]^], then we say that k is soft if A <C 1, and we say that k is collinear to minus if 
77^ ^ 1. Hence, p is collinear to plus and p is collinear to minus in the limit m^/P^ — >• 0. 

2.2 Evaluation of the diagrams 

Now we calculate the double logarithms that arise from the Feynman diagrams that contribute 
to the NLO QCD corrections to the amplitude. The amplitudes for each diagram contain spin and 



color projectors that put the QQ pairs into states of definite spin and color [18]. When the relative 
momentum of the Q and Q in each chamionium is zero, the spin-singlet and spin-triplet projectors 
are given by 

n,{p,p) = -—^^{p + m,), (2.1) 
2v2m 

n,{p,p,X) = + (2.2) 

2v2m 

where £*(A) is the polarization vector for the QQ pair in the spin-triplet state. 

The NLO diagrams that contain the double logarithms in Q^/ml are shown in Fig. |l|. Because 
the process e^e^ J/y+tIc does not satisfy quark-helicity conservation, its amplitude must con- 
tain a numerator factor m^, which produces a helicity flip. Hence, the amplitude is suppressed by a 
factor of nic/ Q relative to a helicity-conserving amplitude. The factor can come from the numer- 
ators of the quaik propagators or from the numerators of the spin projectors. We must also retain a 
nonzero quark mass in denominators because logarithmic collineai^ and endpoint divergences that 
appear in the calculation are sensitive to that mass. 
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(f) 



(b) 



(c) 



(h) 



(d) 





(e) 



(j) 



Figure 1: One-loop diagrams that produce double logarithms of Q^/in^. The upper QQ pair corresponds to 
the //i/, and the lower QQ pair corresponds to the rjc. 



A straightforward calculation of the diagram in Fig. |T](a) gives 



+ 



+ 



Cf - 2<^'4 ) X 



1 



(2;r)^ (A:2 + /£) [{k + p)2 - m2 + /g] [{k - pY - + ie] 
1 



(271)^' [(^ + p + PY + ie] [{k + 2p + pf - + ie] [{k + pY - + ie] 
d''k 1 



(271)^ {k^ + /e) [(fc + pY -nY + ie] [{k + p + pY + ie 



+ • 



(2.3) 



where we retain only the terms that contain the double logarithms in /m^^. Here, Cf = (A'^^ — 
\)/{2Nc), Ca = Nc, d = 4 — 2e, and is the LO (order-a^) contribution to the amphtude 

i^if ^QQiCS,) + QQiCSo)]: 



(2.4) 



where jU. is the polarization of the virtual photon. We use the nonrelativistic normalization for the 
spinors. e^^'^P is the totally antisymmetric tensor in 4 dimensions, for which we use the convention 
£oi23 = 1- We evaluate the first integral in Eq. (2.3) using dimensional regularization to control the 
soft divergence. We find that 



d''k 



1 



(271)^ + ie)[{k + pY -ml + ie] [{k - pY - + ie] 



£lR 



■log(m^/M^) 



log(m2/Q2) + _log2(^2/g2)^. 



(2.5) 



where we retain only the terms that are singular in e and the double logarithms. 
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We have canied out a detailed analysis of the first integral in Eq. ( pj| ) that makes use contour 
integration for the integration over ko. That analysis shows that double logarithm in this integral 
comes from the region in which the gluon with momentum k is simultaneously soft and collineai\ 
That is, it is a Sudakov double logarithm. The second integral in Eq. is identical to the first 
integral if we change the loop momentum to k + p + p. Hence, the second integral gives a Sudakov 
double logarithm that comes from the region in which the gluon with the momentum k + p + p is 
simultaneously soft and coUinear. The last integral in Eq. (23) yields 

f d k 1 ^ri2/2/ 

^ 1 (271)4 (^2 + [(^ + _ ,„2 + [{k + p + p)2 + ie] = ^ °^ ^"''^ 

(2.6) 

where again we show only the double-logarithmic contributions. (In this integral, there is no sin- 
gularity in £.) A detailed analysis, which makes use of contour integration for the ko integration, 
shows that the double logarithm in this integral comes from the momentum region in which the 
gluons carry almost all of the collinear-to-plus and coUinear-to-minus momenta from the spectator- 
quark line to the active-quark line. We call a contribution from this momentum region an endpoint 
contribution. An important observation is that can be made to look similar to the integral that 
gives the Sudakov double logarithm by changing the loop momentum to the spectator-quark mo- 
mentum i = —k — p. Then, we have 

J (271)4 (^2 _ ^2 + + pY + [(^ _ py + • ( •'7) 

It follows that the endpoint double logarithm arises from the region in which the momentum of 
the spectator quark, I, is soft and collinear. From Eq. (2/7), we see that the soft and collinear 
divergences are regulated by m^. An analysis of Eq. (2/7) also shows that single logarithms of 
Q^/m^ can arise from the region in which the momentum of the spectator quark is soft. 

Carrying out similar analyses of the remaining diagrams that contribute to the amplitude 
Y* — > QQ\Q^\-) + QQii^^o) order a^, we find that the double logarithms in each diagram are 
accounted for by Sudakov double logarithms and endpoint double logarithms. The Sudakov and 
endpoint double logarithms that arise from each diagram are summarized in Table. |l|. Our result for 
the sum of double logarithms in all of the NLO diagrams agrees with the results in Refs. [12, p^ ]. 
We also find that our results in the Feynman gauge for the double logarithm in each diagram agree 
with the results that were obtained in carrying out the calculation of Ref. [p^.^ 

Our detailed calculations are consistent with two general properties of the Sudakov and end- 
point singular regions: (1) the Sudakov double logarithms cancel in the sum over diagrams; (2) the 
endpoint region produces only logarithmic singularities, not power singularities. In the following 
sections, we show how these general properties arise from soft-coUinear approximations that are 
valid in the Sudakov and endpoint regions. 



3. General analysis of the Sudakov double logarithms 

As we have mentioned, Sudakov double logarithms arise from a region in which the momen- 
tum of a gluon is simultaneously soft and collinear. We would like to apply a collinear approxima- 

'We thank Yu Jia and Xiu-Ting Yang for providing us with the logarithmic contributions of the individual NLO 
Feynman diagrams. 
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Diagram 


Endpoint double logaiithm 


Sudakov double logarithm 


(a) 


{Cf — jCa)<S' 


l{Cf - \CA)y 


(b) 


CfS 





(c) 


ICfS 





(d) 


^ CfS 





(e) 


{Cf — \Ca)S 


{Cf - \CA)y 


(f) 





{Cf - \Ca)^ 


(g) 





-{Cf-\CA)y 


(h) 





— {Cf — \Ca)'^ 


(i) 





— {Cf — \Ca)'9' 


0) 





—{Cf — \cA)y 



Table 1: Endpoint and Sudakov double logarithms in each diagram, in units of /jz^Io x {—iOLsTtQ^) /2. 



tion to such gluons. Consider, for example. Figs. |T](e) and (g), in which a gluon with momentum 
k that is collinear to minus is emitted from a quark line with momentum p. Then, the quark-gluon 
vertex and the two propagator numerators surrounding it can be written as 



{p + mc)Y^'{p + fc + m,)=2{p^' +k^'){p+m,)-fcy^m. 



(3.1) 



where /i is the polarization index of the gluon, we have used fcp oc pp and pp = ml, and we have 
dropped terms of order m^. In the coUinear-to-minus approximation, one retains only the first of the 



two terms on the right side of Eq. (3.1). In the case of nonzero quark masses, this approximation is 
not valid in general. However, if k is soft in comparison to p, as well as collinear, then we can drop 



the second temi on the right side of Eq. (3.1), and the standard collinear approximation holds. [We 



can also drop k in the first term on the right side of Eq. (3.1).] Since the cuixent in Eq. (3.1) now 
lies in the minus light-cone direction, up to terms of order ml, we can make a collinear-to-minus 



approximation in the gluon propagator [ ]19 
tensor 



guv 



|, by making the replacement in the polarization 



k- p — ie 



(3.2) 



where the index v corresponds to the attachment of the gluon to the quark line that is collinear to 
minus and the sign of is is fixed by the sign in the original Feynman diagram. This approximation 
is valid unless the jj. attachment of the gluon is to a line that is also collinear to minus. Hence, 
the approximation always holds for the diagrams that produce Sudakov logarithms because the 
invariant in the logarithm can appear only if the soft-collinear gluon connects a line carrying 



momentum p with a line carrying momentum p. The replacement (3.2) can also be regarded as a 
soft approximation ||22, 23] to the n attachment of the gluon. However, the collinear approximation 
can be more useful in applications other than the present one because its form is independent of the 
direction of the momentum p, while the form of the soft approximation is not. 



For the diagram of Fig. [l](e) we can apply the soft-collinear approximation (3.2), where the 
index jj. corresponds to the connection of the collinear-to-minus gluon to the active-quark Mne that 
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canies momentum p: 

1 in^ 1 

u{p)fi—, -^uip) , ■ J. f - = -u{p) , . , (3.3) 

where, in the last equahty, we have apphed the graphical Ward identity (Feynman identity). Sim- 
ilarly, for the diagram of Fig. |](g) we can apply the soft-colUnear approximation (^^, where the 
index pL corresponds to the connection of the coUinear-to-minus gluon to the spectator-quark line 
that canies momentum p: 

1 1 /^D^ 

— ^-^ -y^vip) ^ — ^-^ . v{p) = + , . v{p). (3.4) 

-p+Ji-ruc + ie' ^ ' -p+Ji-Mc + iek-p-iE ^ ' k-p-ie ^ ' 

Because the outgoing QQ pairs are in color-singlet states, these two contributions have the same 
color factor and cancel. This type of cancellation extends to all diagrams involving a gluon with 
soft-coUinear-to-minus momentum. Similar cancellations occur for diagrams involving a gluon 
with soft-collinear-to-plus momentum. Therefore, in the sum of all diagrams, Sudakov double 
logarithms cancel. 



4. General analysis of the endpoint region 

As we have mentioned, the endpoint double logarithms of /m^. arise from the region of 
loop integration in which the momentum i of the internal spectator-quark line is simultaneously 
soft and collinear. Hence, we need to consider only diagrams that can produce such a momentum 
configuration. These diagrams are shown in Figs. |](a)-(f). (A diagram and its charge conjugate 
give equal contributions to the amplitude. The charge-conjugate diagrams are not shown in Fig. [I|) 
The endpoint double logarithms arise from contributions in which there is a numerator factor nic, 
which produces the helicity flip, and in which integrals diverge logarithmically in the limit nic — 0. 
In general, integrals can also diverge as inverse powers of in the limit nic — ^ 0, but, as we shall 
see, such contributions vanish when the numerator trace is taken. 

In diagrams (a), (e), and (f), the momenta of the propagators on the active-quark lines contain 
both p and p. Since p- p ^ P'^ ^ Q^, we can ignore I in the denominators of those propagators. In 
the limit nic — ^ 0, the two gluon-propagator denominators and spectator-quark-propagator denomi- 
nator produce factors 1 / + 2p • £) , — 2p-l), and respectively, where we have dropped 
the ni^ terms in the propagator denominators. Hence, in order to obtain a logarithmically divergent 
soft power count (A^"^) and logarithmically divergent collinear power counts [(t]*)^"*], we must 
drop all numerator factors of i. This implies that the helicity flip comes from the factor nic in the 
numerator of the spectator-quark propagator. 

In diagram (b), the momentum of the outermost active-quaik propagator contains the momen- 
tum p, but not the momentum p. Hence, in the limit nic — ^ 0, the denominator of this propagator 
produces a factor \ / {l^ + 4p ■ I). Then, all of the propagator denominators taken together pro- 
duce a linearly divergent soft power count and a linearly divergent coUinear-to-plus power count. 
However, it is easy to see that numerator factors reduce both of these power counts to logarithmic 
ones. First, we rewrite the numerator factors that are associated with the outermost gluon and the 
spin projector for the QQ pair as y^{p — nic)^*y^ = 2nic^*, where we have used the fact that 
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p-e* =0. Now, because the numerator power of m^. is in this factor, the numerator of the spectator- 
quark propagator must contribute a factor /. Furthemiore, the numerator vanishes, up to terms of 
order m^, if ^ is proportional to p because there are then two factors of that are separated by 7 
matrices with which they anticommute. Hence, in the trace over the 7 matrices, £ must appear in 
the combination i ■ p. This numerator factor reduces both the soft and collinear power counts to 
logarithmic ones. 

In the case of diagram (c), the denominator of the outermost active-quark propagator con- 
tributes a factor —4p-i). Hence, the soft and coUinear-to-minus power counts from the prop- 
agator denominators are both linearly divergent. Now, we rewrite the numerator factors that are as- 
sociated with the outermost gluon and the spin projector for the Tjc QQ pair as Y^{—p — nic)Y5Y^ = 
{—2p + 4mc)Y5. If the numerator factor comes from the spectator-quark propagator, then there 
must be a factor / from the outemiost active-quark propagator, or else there are two adjacent fac- 
tors of p and the numerator vanishes, up to teniis of order m^. If the numerator factor nic does not 
come from the spectator-quark propagator, then that propagator must yield a numerator factor /. 
In both cases, the numerator factor vanishes, up to terms of order m^, if i is proportional to p be- 
cause, in that case, there ai^e two factors of p that are either adjacent or are sepai^ated by 7 matrices 
with which they anticommute. Therefore, we conclude that the trace contains a factor i • p, which 
reduces both the soft and collinear-to-minus power counts to logarithmic ones. 

In the case of diagram (d), the denominators of the active-quark propagators produce factors 
I / + 4p ■ tj and \/{f' —Ap- tj. Hence, the propagator denominators, taken together, produce a 
quadratically divergent soft power count and linearly divergent collinear-to-plus and collinear-to- 
minus power counts. One can apply the arguments that were used for the numerators of diagrams 
(b) and (c) separately to each of the gluons in diagram (d), with the conclusion that the numerator 
contains two factors of / and that the numerator vanishes, up to terms of order m^, if I is proportional 
to p or to p. Hence, the trace contains a factor £■ p£- p or a. factor In either case, the numerator 
factor reduces the soft and collinear power counts to logarithmic ones. 

In the case of one-loop collections to helicity-consei^ving chamionium-production processes, 
in which there is no numerator factor m^, the previous arguments show that the contribution from 
the region in which the spectator quark carries a soft-collineai^ momentum vanishes, implying that 
there are no endpoint double logarithms. This general conclusion is confirmed in explicit calcula- 
tions [M|24|, ^]. 



5. Summary 

In this work we have investigated the origin of the double logarithms of Q^/ml that appear in 
the NLO QCD corrections to the process e+e^ J / y + Tic- We find that the double logarithms 
in each diagram are accounted for by Sudakov double logarithms and endpoint double logarithms. 
The Sudakov double logarithms cancel in the sum of all diagrams. We have given a general ar- 
gument for this cancellation that is based on the soft-coUinear approximation and graphical Ward 
identities. We have reinterpreted the region of a loop integration that gives rise to an endpoint 
double logarithm as a leading region in which the momentum of the spectator-quark line is both 
soft and collineai\ Under this reinterpretation, we would also expect single logarithms of Q^/ml 
to arise from a leading region in which the momentum of the spectator-quark line is soft. This 
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reinteipretation may be useful in establishing an all-orders factorization theorem for helicity-flip 
quarkonium production. Such a factorization theorem might allow one to resum logarithms of 
Q^/ml to all orders in for helicity-flip processes. We have also given a power-counting analysis 
of the one-loop endpoint contributions, which shows that they can arise only in the presence of a 
helicity flip and that the loop integration can produce logarithms, but not inverse powers, of the 
heavy-quark mass Mc. 
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